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We quantize tlie exterior of splierically symmetric vacuum space-times using a midi-superspace 
reduction witliin tlie Asiitekar new variables. Tlirougli a partial gauge fixing we eliminate the 
difTeomorphism constraint and are left with a Hamiltonian constraint that is first class. We complete 
the quantization in the loop representation. We also use the model to discuss the issues that will 
arise in more general contexts in the "uniform discretization" approach to the dynamics. 



I. INTRODUCTION 

00 ' 

' Loop quantum gravity has emerged in recent years as a significant candidate for a theory of quantum gravity. See 
. [l[ for recent reviews. The theory has a mathematically rigorous basis for its quantum kinematics [2], which has also 
^sj ' been proven to be unique up to certain assumptions Q . The problem of the dynamics of the full theory has remained 
sj^' unsettled. The origin of the difficulties can be traced back to the kinematical space of states. In this space, infinitesimal 
' diffeomorphisms are not implemented as operators. The issue is further compounded by the fact that one cannot 
I check the consistency of the quantum constraint algebra, since the commutator of two Hamiltonians is proportional to 
■ an infinitesimal diffeomorphism. Attempts to circumvent this problem, like representing the Hamiltonian constraint 
00 ' as an operator on states invariant under diffeomorphisms, are faced with other difficulties (see [4] for a more thorough 
discussion) . 

' These problems have led several researchers to consider alternatives to the usual Dirac approach to the problem 
of the dynamics. One of the alternatives is the "master constraint" project of Thiemann and collaborators [5] which 

• has similarities to an earlier proposal by Klauder Q. Others consider the covariant "spin foam" approach as an 
I ] alternative, since one may bypass the construction of the canonical algebra entirely, at least in some settings. Another 

• approach is the one we have presented in recent papers called "uniform discretizations" [7, 8]. 

' The spherically symmetric case has also specific problems in the traditional approach to loop quantum gravity. 
C ■ i It has not been possible up to now to find a particularization of the construction of Thiemann for the Hamiltonian 
' constraint to the spherically symmetric case that has the appropriate algebra of constraints on the diffeomorphism 
Q , invariant space of states [9] . 

O^' In spite of these difficulties, the loop approach has been successful in the context of homogeneous cosmologies, giving 
^ '] rise to "loop quantum cosmology" ■ The reason why the approach works in this context is that there is no issue 
with the algebra of constraints: the diffeomorphism constraint is gauge fixed by the use of manifestly homogeneous 
^ variables and there is only one Hamiltonian constraint. We would like to show that a similar situation emerges in 
the case of spherically symmetric space-times. We will show that one can fix the diffeomorphism gauge in such a way 
that one is left with only a set of Abelian Hamiltonian constraints. 

Since this is our first approach to this problem, we will not explore the more interesting possibilities of this model, 
for instance, what happens to the singularity. We will restrict our attention to the exterior of the horizon. We will 
also fix a gauge to eliminate the diffeomorphism constraint for the sake of simplicity. The resulting Hamiltonian 
constraint has nevertheless a non-trivial first class algebra (with structure functions) in the continuum, but we will 
show it can be Abelianized. We will show that the loop quantization can be completed within the traditional Dirac 
quantization approach. In particular we will construct the kinematical and physical Hilbert spaces and we will show 
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the quantization to be equivalent to the one carried out in terms of traditional variables by Kuchaf [llj. We will 
discuss the implications for the elimination of the singularity, though we will not work out the details in this paper. 

We would also like to discuss what lessons one can get from this model for the "uniform discretization" approach to 
the dynamics of quantum gravity. Since the uniform discretization approach is equivalent to the Dirac quantization 
procedure if one has Abelian constraints, it does not add anything new to the quantization of this model if one chooses 
to Abelianize the constraints as we do. We nevertheless would like to discuss what would happen if one had chosen 
not to Abelianize the constraints. The uniform discretization procedure can be applied, but the treatment of the 
model becomes more involved, and the degree of complexity increases with how much one chooses to ignore about 
the particular details of the model. 

The organization of this paper is as follows. In section II we set up the classical variables for spherically symmetric 
space-times in the Ashtekar formulation; in section HI we discuss the loop quantization using the traditional Dirac 
procedure and in section IV we discuss the use of uniform discretizations. We end with a discussion. 
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II. SPHERICALLY SYMMETRIC SPACE-TIMES 



We will use the Ashtekar new variables to describe the spherically symmetric space-times. Previous work on 
this subject was done by Bengtsson [l^, Thiemann and Kastrup [l3| which also Abelianize the constraints with 
the traditional (complex) Ashtekar variables and with the more modern real connection variables by Bojowald and 
Kastrup [U, and Bojowald and Swiderski Q so we will not repeat the full construction of spherically symmetric triads 
and connections (the latter were also discussed in the context of spherically symmetric Yang-Mills theory by Cordero 
and Teitelboin i |15l| ). Also, in the context of geometrodynamical variables, gauge fixings in spherical symmetry were 
considered by [la]. 

One assumes that the topology of the spatial manifold is of the form E = i?+ x S^. We will choose a radial 
coordinate x and study the theory in the range [0,cxd]. We will later assume that there is a horizon at x = 0, with 
appropriate boundary conditions as we discuss below. 

The invariant connection can be written as, 

A = A^{x)A3dx + {Ai{x)Ai + A2{x)A2) dO + ((Ai(x)A2 - A2{x)Ai) sin 6 + A3 cos 6) dip, (1) 

where Ax, Ai and A2 are real arbitrary functions on the A/ are generators of su{2), for instance A/ = — icr//2 = r/ 
where aj are the Pauli matrices or rigid rotations thereof. The invariant triad takes the form, 

o /-) Q 

E = E-{x)A3sme— + {E'{x)Ai + E\x)A2) sin0— + {E^{x)A2 - E^x)Ai) — , (2) 

where again, E^ , E^ and E^ are functions on One has the following canonical Poisson brackets for the symmetry 
reduced variables, 

{Axix),E^x')} = 2jG6{x-x'), (3) 
{A,ix),E\x')} = jG6ix-x'), (4) 
{A2{x),E^{x')} = jG6{x-x'), (5) 

and all the other brackets vanish; G is Newton's constant and 7 is the Immirzi parameter. 

Under coordinate transformations of the x coordinate x'{x), A^ transforms like a scalar density A'^{x') = 
{dx/dx')Ax(x) whereas E^ is a scalar. A^ and A^ are scalars and E^^ and E^ are scalar densities. 

There are three constraints. The first one is a Gauss law, 

G{\) = j dxX {{E'^y + 2AiE^ - 2A2E^) , (6) 
that generates U{\) gauge transformations on the line. The second one is the remnant of the diffeomorphism constraint, 

D[Nx\ - j dxNx (2(Ai)'i?i + 2{A2)'E^ - Ax{E-)') , (7) 
and finally the Hamiltonian constraint, 

H[N] = {2G)-^jdxN{x){\E''\[{E^f + {E'^f]y^''^x{2E'-'{E^A'2-E'^A\)+2AxE^{AiE^+A2E'^) (8) 
+ {Al +AI-1) {{E^f + {E^f) - (1 + 7') {2KxE^ {K,E^ + K2E^) + {{K^f + {K^f) {{E^f + [E^f)) , 
where K.^, Ki and K2 are the independent components of the spherically symmetric curvature, 

K = KxA^dx + (i^iAi + K2A2)d6 + {K1A2 - K2A1) sinOdip, (9) 

which can be written in terms of the canonical variables A's and £"s. 

It simplifies things if one introduces "polar" type coordinates in the directions 1, 2. To do this we define, 

Al = A^cos/3, (10) 
A2 = A^s\a(3. (11) 

We introduce a canonical transformation from the variables Ai, A2, E^ , E^ to the variables A^,(3, P'^ , through the 
type III generating function F = E^A^ cos (3 + E'^A^p sin /3. This leads to the above relation between Ai^2 and A^p and 
/3 and defines the conjugate momenta, 

P'^ = 2E\os^ + 2E'^smP, (12) 
P>^ = -2E^A^sm(3 + 2E'^A^cos(3, (13) 
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and introducing a new angular variable a via = E'^ cos (a + /3) and E"^ = E'^ sin (« + /?), we have that P"^ = 
2i5'^ cos a and = 2EVA^ sin a, and we have rescaled P'^ and P^ by a factor of two so one has the canonical 
Poisson bracket relations, 

{A^{x),P^{x')} = 2-iG5{x-x'), (14) 

{P{x),P^{x')] = ^^GSix-x'), (15) 

with the symplectic structure in Ax, E^ unchanged from before. 

In terms of these variables the Gauss law and diffeomorphism constraint simplify, 

G = P'3 + (£;^)', (16) 

D = Pl^ fi' + Pf A'^ - {E'=)' Ax. (17) 

To write the Hamiltonian constraint it turns out that a further change is desirable towards variables more closely 
connected with the geometry. Let us start by identifying the metric constructed from the densitized triads we have. 
To do this we first write the determinant of the metric deig = \E^\{E^)'^ svc? 9. The metric components are, 

fc = ^|^, 9eo = \E% g^^^\E^\s\n^e. (18) 

We consider generators of su{2) rotated with respect to the Cartesian basis. Given tan(a + /3) = —E'^/E^, we define, 

A| = Aicos(a + /3)+A2sin(a + /3), (19) 

A| = -Aisin(a + /3) + A2Cos(a + /3), (20) 

so the (undensitized) co-triad can be written as, 

e = e^Aada; + elK%d9 + e^Kld^, (21) 

where, 

e3 = - (22) 

el = VWl (23) 

e% = y/\E^\sine. (24) 

From the co-triad we can compute the spin connection, 

r3=-(a + /3)', r^ = ^, r^ = -Msin0, Tl=cose, (25) 

e-x ^x 

and then from jK^ = A\^ — we can compute the curvature components 

^Kx^Ax + [a + P)', (26) 

7ifi = Ai-F^sin(a-f/3), (27) 

7/^2 = ^2 + cos(a fi) , (28) 

where we define F;^ as, 



2EV 



r.--^ = -^, (29) 



and also F^, = Fi^. 

As was noted by Bojowald and Swiderski [9], in order to introduce a representation in terms of holonomies one 
requires certain falloff conditions, in particular A^p — > as x — s- oo. This unfortunately is not true. To see it, note 
that F;^ is a ratio of two densities and therefore a scalar. But asymptotically, E^ ~ x"^ and therefore (E^)' ~ 2x 
and E''^ ~ x where M is the mass of the classical solution considered and therefore F;^ ~ 1, which implies that A^p 
will not tend to zero asymptotically since K^p asymptotically and we have the relation {A^)'^ — (F^)^ + {jK^)'^ 
{K^p = Kl + K2 ) . For further discussion of asymptotic properties see Kuchai^ (111] . 
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To construct a connection with a good asymptotic behavior we consider the following canonical transformation, 

A^^ = 2 cos a A^, (30) 

p^rj = a + /3, (31) 

with the following type II generating function 

F ^ P'>{a + P) + 2E'^A^cosa, (32) 

which recovers the above transformation and introduces the canonical momenta P'' and , 

P>^ = P'', Pf = 2Efcosa, (33) 

which leads us to consider the following canonical variables A^p, E"^ ,7], P^ , in terms of which the Gauss law and 
diffeomorphism constraint take a simple form, 

G = P" + (£:^)' (34) 
D = P'^T]' + E'^A'^- (E'^YA^. (35) 

Before continuing it is worthwhile noting that E^^E"^ are "metric" variables, in the sense that they are invariant 
under the transformations generated by Gauss' law. A^p and A^ + rj' are also invariant and the latter is proportional 
to Kx- In order to clarify the meaning of A^p notice first that, 

yl^A^ = A^cos/3A2 -^y,sin/3Ai, (36) 
= ^iA2 - A2hi = r^A| + iK^Kl, (37) 

and then note that A|; and A*^ are orthogonal, i.e., Tr(A|;A^) = and also that Tr(A^A^) = cos a. Therefore 
A^p cos a — J Kip and it follows that A^p — 2^Kp. This automatically implies that the connection has the right falloff 
condition since Kp, 0. We also have that, 

-fK, = A,-T,^ A., + [a + p)' = Ax + ri\ (38) 

which confirms what had been noted, that A^ + 77' is a gauge invariant combination. 

It is convenient to relate these variables of the Ashtekar formalism with the canonical variables used for the spherical 
case by Kuchaf . The latter are given by, 

ds^ = k{xfdx^ + R{xfdn^, (39) 
with A > and R > and the latter is the curvature of the 2-spheres. We then have, 

K,^ = -AN-' (a - (iV^A)') = -AK, + 2^i^^ (41) 
Kgg ^ -N-'R{R- R'N'') ^ ~RK^ (42) 
One therefore has a canonical transformation between the pairs {Ap, E^ , A^ + r/, E^) and (Pa, A, Pr, R) where, 

Pa - Vm^, (44) 
27 

R = (45) 
-A^ + r( E'P A. 



A - (43) 



P« = (46) 



7 .JW\ 27 

The Hamiltonian constraint can be obtained from equation ([5]), and takes the form 

E^ A^A^.JW\ , {{E'^yf 



H 



2^/\Ey\ 272 8y|^^2 %,J\fF\Ev 

^\{E^y{E^y _ A^y/WW , ^/W\iEn" .47. 

2{E-fi)^ 272 2E'P ' ^ ' 
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We will now fix the spatial difFeomorphism gauge freedom. This will simplify calculations considerably. We would 
like to have two things: a) that the Gauss law remains as a constraint in order to have usual loop representation 
techniques for quantization and b) that the expression for the spatial volume retains a simple form since it plays an 
important role geometrically. We choose a gauge that is suitable for the region exterior to a horizon, = {x + a)'^ 
where a is a positive constant given by Horizon = cl^, which is equivalent to R = x + a and x — corresponds to 
the horizon. This gauge choice commutes with Gauss' law, which therefore remains a first class constraint. Since we 
do not know the value of E^ at the horizon, a is really a dynamical variable, its conjugate momentum is related to 
x^Q. In that gauge the diffeomorphism constraint takes the form, 

E^A'^ - {E^'iA, + 77') - E'^A'^ - 2{x + a){A, + r^') = 0, (48) 

and it can be explicitly solved for Ax + rj' — E'^A'^/ (2 (2; + a)). This also determines strongly the value of the canonical 
momentum of the variable a. 

After solving for that variable, the Hamiltonian constraint takes the form, 

Ef (Al{x + a)\' Ef 3{x + a) , f I \' 

^--(^-r^(^V^j -2(^ + ^ + ^-+«) U) 

The constraint only depends on the canonical pair Aip{x), E^ (x) and a. The constraints remaining after eliminating 
the diffeomorphism constraint are first class, but have an algebra with structure functions, 

My)„^,\'u ^ My) 



{H{x),H{y)} = ^-^H{y)J 5{x ~ y) ~ -^H(y)5,x{x ~ y), (50) 

{G{x),H{y)} = 0, (51) 
{G(x),G(y)} = 0. (52) 

This reduced system after eliminating the diffeomorphism constraint is the one we will use as starting point for the 
quantization. 

Our gauge choice is suitable for describing the exterior of a black hole. The initial data for the system is given on a 
spatial surface extending from the Schwarzschild radius to infinity and it determines entirely the spacetime bounded 
by the horizon and null infinity, since no entering data is given on these two surfaces. The boundary conditions on 
the horizon correspond to having the gxx = {E"^)"^ /{^ + a)^ component be singular, l/E'^\x=o = and A^|a;=o — 
(as discussed by Bojowald and Swiderski^y'l ) which corresponds to the isolated horizon boundary condition For 
the falloff conditions at infinity we have , 

E"^ = x + M{t) + 0{x-'), (53) 
A^ = 0(:e-(i+^)), (54) 

where M{t) is a function of t and the above expressions hold in the case that the value of a; M{t). One also has 
that the shift vanishes asymptotically — 0{x~'^) and the lapse N = N + 0{x^'^). 

In order for the variations of the dynamical variables to preserve the falloff conditions one has to add boundary 
terms to the action. Given, 

S{A^,E'^,N,a) ^ lim j dt j^^ dx {e^ A^ - N H{x)^ . (55) 

The problematic term is the one stemming from the variation of E"^ in the term (x + a)^(l/i?'^)' in the Hamiltonian. 
Noting that SE"^ — 5M+ and that in problematic term one has, 

dxN{x + af {{E'^y^SE^y) . (56) 

This term requires that one integrate by parts. This produces an extra term of the form N{xjf-)SMjf- and it would 
require that the lapse vanish at infinity. This led Kuchaf to propose adding a term at infinity of the form — J dtM^f+ 
with T+ a new dynamical variable such that its variation yields Af_|_ = and variation with respect to E'^ is now well 
defined and yields f+ = A^+. The resulting action is, 

S{A^,E'^,N,a,T)^ lim J dt J^^ dx (^E'^A^ ~ NH{x)^ - J dtM+T+. (57) 

Through the partial gauge fixing we are now left with a model that is still a midi-superspace but that has only one 
constraint, the Hamiltonian, with a non-trivial first class algebra. 
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III. QUANTIZATION 

Our objective is to use the spherically symmetric model to explore the properties of the loop quantization approach 
in a case with field theoretic variables. Although the model in the end has a finite number of degrees of freedom, this 
is achieved in a non-trivial way through the imposition of the Hamiltonian constraint, all the time treating the model 
as a field theory. 

We will for simplicity concentrate on the exterior of the black hole spacetime. We will therefore not concern 
ourselves in this first approach on issues like how the use of loop variables may eliminate the singularity, etc. The 
issues we are interested in probing are: the role of the Hamiltonian constraint in generating the evolution, seeing if 
the theory is coordinate invariant after quantization and discretization and how the loop treatment compares with 
the continuum quantization and seeing if one can avoid dealing with anomalies in the quantum constraint algebra. 



A. Traditional canonical quantization 

In order to compare with the other quantization approaches, we briefly review the traditional canonical quantization 
of this model, first discussed by Kuchai^ 11] (see this paper for earlier reference on this subject). In the region x > 
exterior to the horizon, one can fix a gauge globally. For instance = 2^K^p = 0. The constraints A^p = and 
H(x) = are second class, since their Poisson bracket is non- vanishing (not even weakly). To quantize one has 
to impose these constraints strongly. In the action the only remaining contribution after this gauge fixing is an 
asymptotic one, 

S = j dtfM+. (58) 

From here one gets canonically conjugate variables r, P'^ = — Af+. t is the proper time that determines the position 
of the spatial hypersurfaces of vanishing extrinsic curvature (usual Schwarzschild slicings). For these surfaces one has 
A^p = and the quantity E"^ is the solution of. 



{ x + g)^^; 3{x + a) _ 

{Evy ^ 2EV 2{x + ay 



H = = - ' ^' V ^ -^y-^-^^' ^ ^ (59) 



which yields as solution 



E^ = + = A_ (60) 

and one can identify by expanding asymptotically that M — a/2 and recalling that gxx = {E'^Y jE^ one obtains the 
usual form of the Schwarzschild metric with M the mass of the space-time. 

The quantization is straightforward, since the only remaining canonical variables are M and r. These variables 



have no dynamics. One has A^ — 0, E'^ = (a; -I- 2M)/ y 1 ^ 2M /{x + 2M) and one can for instance introduce an 

eigenbasis of the mass operator, M(j)(m) = •m<j){m). 

The operators associated to the other components of the metric can be determined easily. The lapse is determined 
by the equation fixing the preservation in time of the gauge condition, A^ — {A^, J N(x)H{x)dx} = with solution 
N = y^l — a/{x + a). The shift is determined by noting that if A^p = then — and therefore equation ([42]) 
implies, given that R — 0, that = One then has goo = N'^ and this leads to the usual Schwarzschild solution: 
300 = (1 - 2M/R). 

Different choices of gauge lead to quantum theories in different coordinate systems. In the gauge chosen above 
A^ = 0. For an more general gauge with arbitrary A^p one has that 



E'P ^{x + 2M) 



1 - 



2M A, 



2M 472 



(61) 



For instance in the gauge A^p = ^k,x/{x -f 2M)2, with k a constant, and which satisfies the boundary condition at 
the horizon and at infinity, one has that, 

X + 2M , , 

E"^ = , (62) 

/1 _ 2M I KX 

X+2M i(x+2M) 
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and the lapse and shift are both non- vanishing, 



= y-^T^ + Aix + 2Mr V^^T^J ' ^^^^ 

NA 

= (64) 

In the end all the quantizations are equivalent since in essence we are dealing with a mechanical system and the 
choices of coordinates are just different choices of quantities that one computes for the mechanical system, and one 
has only two observables, the mass and the proper time at infinity. 



B. Loop quantization 

In the loop representation the fundamental operators are associated with holonomies. This requires re-casting 
the theory of interest in terms of such variables. This usually involves taking limits of holonomies around loops of 
vanishing area, for instance, to represent the curvature. Moreover, the operators we wish to consider need to be 
regularized. This is usually achieved by discretizing the theory on a lattice. Therefore the usual loop treatment of 
theories is a natural framework in which to apply the "uniform discretization" technique. 

We will proceed as follows: 

1) We will identify the space of states of the loop representation for the spherically symmetric case and the operators 
that are well defined in this space. 

2) We will particularize the space of states to the gauge we are considering in this paper, where one is left with 
only the Hamiltonian constraint. 

3) We will note that one can Abelianize the constraint and, upon discretization, the resulting discrete theory can 
be treated using the traditional Dirac quantization procedure (in the case of Abelian constraints it is known to be 
equivalent to the uniform discretization procedure). 



C. Tlie space of states in the loop representation for spherical symmetry 

The quantization in the loop representation is based on cylindrical functions that depend on the connection through 
"open holonomies" . The latter are associated with graphs on the spatial manifold (in this case the one dimensional 
line) composed of a set of edges without intersections g — [J - Ci where g is the graph and are the edges. The edges' 
vertices form a set V{g) composed by all the endpoints of e^. The Ay, correspond to directions transverse to the 
radial one and are represented in the loop representation via "point holonomies" exjp{ifj,yAip(x)) — e:xjp(2ij^yK^) at 
each vertex v with /i G K. as is usually done for scalar fields [l8|. The variable 77 is an angle and the corresponding 
point holonomy is given by exp (invrj{v)) with riy G Z. Before imposing gauge invariance under the transformations 
generated by Gauss' law the spin network states in a basis of spherically symmetric connections are, 

^g,fe,n,M(^) = II*^^P ( ^'^W ^^^^) n exp(i/i.„7i4:^(w))exp(m„?7('y)) (65) 

where fee G ^ is the multiplicity of the loop. If we recall that the Gauss law is + [E^)' — and that A^ and 
Ax + 77' are gauge invariant, one has that the gauge invariant spin networks for the spherically symmetric case are, 

^9,kM) = n {\^^ I + '?') n (2*M.7^^(«)) • (66) 

In the full three dimensional case the Hamiltonian constraint is written in terms of holonomies along small closed 
paths. One can do something similar here for the holonomies along the radial direction, 

exp (i [ ' A^di^ = 1 + ieAx{v) + 0{(^), (67) 



where e — Vf ~ Vo 'is the distance between two vertices connected by an edge. For the directions where we have 
homogeneity we cannot use this argument since we do not have edges that we can make "small" . In that case we can 
seek an alternative: to expand the "point holonomies" in the limit in which p is small, 

exp {i-ipK^v)) = 1 + i-ipK^ + 0(p2) . . . (68) 
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It should be noted that quantum mechanically the limits e ^ and p ^ end up not being well defined and the 
resulting quantum theories therefore approximate the classical theory well in regions where the relevant variables in 
the exponent are small in order to make the above expansions accurate. 

At this point it is worthwhile mentioning that in the full 3D case, in the loop representation, one represents the 
Hamiltonian constraint on the space of diffcomorphism invariant states. Here we will not proceed in such a way. 
To begin with, since we already fixed gauge for the diffeomorphism constraint, there is no sense in representing the 
Hamiltonian in a space of diffeomorphism invariant states. Moreover, if one were to try to mimic the construction 
that is done in the 3D case (without fixing the gauge for diffeomorphisms) , in the spherical case one would encounter 
and additional difficulty: since there is no notion of "planar vertices" the constraints will not close the appropriate 
algebra even if represented on diffeomorphism invariant states, at least with the usual treatment, as discussed in 
Therefore even if one started with diffcomorphism invariant states the theory would not be invariant under space-time 
diffeomorphisms . 



D. Loop representation of the partially gauge fixed theory 

We will work with a space where the spatial coordinates have been gauge fixed in the way we described in the 
classical section, and therefore the diffeomorphism constraint is not present. The relevant canonical variables are 
Aip, E'f (we will usually substitute A^p for since they are simply related) and the asymptotic variables M — a/2, r. 

We will work in a one dimensional lattice L with points 0, . . . x^r in with spacing x^+i — x„ — e(n), and we 
include the possibility of unequal spacing since we will see it may be useful in certain cases. Strictly speaking, for 
the example we are considering, there is no need to discretize and one could work directly in the continuum. We will 
work in a discretized lattice to be able to make better contact later with the "uniform discretization" approach and 
also for comparison with other situations. We will see that in the end we can correctly remove the lattice regulator 
and produce a continuum theory. 

The Hilbert space in the bulk is given by, 

H ^ L"^ ((g)Ari?Bohr, ^Af^Mo) (69) 

where i?Bohr is the Bohr compactification of the real line and /ig is the measure of integration we discuss below. In 
this space one can introduce a basis, 

TgAKv]= n exp{ifi,jK^iv)) (70) 

veV(g) 

with fiy G g the graph, V{g) vertices of 17 in L and the notation /2 denotes (/x„j, . . . ,/it,p) with p the number of 
vertices in the graph. In the gauge fixed case we are considering a graph is just a collection of vertices. We will now 
introduce a notation more adapted to the uniform discretization setting by labeling the points in L with an index 
m,,n, . . .. We would then have K^^n = K^{xn)- We have also rescaled in such a way that {ivTi^.m, £"^1 = GS„i.n 
with G Newton's constant and we choose units where h = c = 1 and therefore G = ^pianck- One then has the quantum 
representation, acting on states ^'[ifi^^m] = {K^^ml"^), 

= ^^^Planck^r; J ('''1) 

and therefore, 

^ni^g,p. ~ ^ ^ A^m7^Planck'^m,n(u)2^g,/Ii (72) 

where the (5 is a Kronecker delta and n{v) is the position on the lattice L of the vertex v of the spin network. 
Given an interval / in i the volume of the corresponding "shell" is given classically by, 

y(/)=47r^ |i?,^J(a:™ + a), (73) 



mel 



and as a quantum operator. 
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We can introduce a basis of loop states \g,il>, 



(75) 



and the Bohr measure guarantees that, 



(76) 

We define the "holonomy" associated with the "transverse" connection at a vertex v. We will use this to construct 
the Hamiltonian. The definition is, 

(77) 



h^{v, p) = exp ( -pA^{v) ) = exp {ip-fK^{v)) , 



one has that, 

/i¥'(w,,p)|.g,M l.9,Mi>i, • • • ,Mi>. +P,--->, (78) 
and in the case of vi not belonging to the initial spin network, the action adds a vertex with fi = p. The co-triad is 

-E'ml5iAt>= XI ^^vl(-l\a.nck.^m,n(v)\9■,^^> ■ (79) 

We will now define the inverse operators that arise in the Hamiltonian. It should be noted that similarly with 
what happ ens in loop qiiantTim cosmology, the inverse operators are bounded. We will start by computing classically 

sgn(^;'^)/^|^'^|. To do it wc note that, 



cos 



(^) {#^,sin (^) } - sin (^) {^,cos (^) } 



GsgnCB^ 2 (jpK^ GsgnCB^ 2 ( IpK^ 



Gsgn(i;;^)7P 



V 2 



which leads to 



Quantum mechanically one has 
sgn(£;^) 



' J-^m. 



9,l^> = 



Ai 



4^ 



\ 2 



Dm. 



7^1 anc 



sm 



-E^l.cos 



|-E^|,sin j 

/ IpK^^rn 



IPK^ 



\9,IJ'>, 



and noting that. 



one finally has, 



^J\^\\g,{i>= J2 ^ 



vevig) 



sgn(g;^) 



(80) 
(81) 
(82) 

(83) 

(84) 
(85) 

(86) 
(87) 



«ey(s) 



The maximum value of this quantity occurs for p.^ = p/2 and is given by 2/(y^£pianck)- For p. ^ p one should 
recover the classical approximation. In that limit \g,il > is an eigenstate with eigenvalue 1 / ( \/7-^Pianck -\/ 1 Mi; I ) and 
satisfies that, 



\Ev\ 



sgn{Ef 

~7W\ 



\9,P') = \9,P). I^»P 



which is the usual relation between the classical variables. In what follows we will either work in the loop representation 
or the connection representation as needed. 
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E. The Hamiltonian constraint, Abelianization and traditional Dirac quantization 



Here we will show that one can Abelianize the Hamiltonian constraint, discretize it and then quantize the discrete 
theory using the traditional Dirac quantization procedure. It is known that for Abelian constraints the uniform 
discretization procedure coincides with the Dirac quantization so we will not work it out explicitly. This quantization 
will provide a baseline against which to compare quantizations in which we do not Abelianize the constraint. The 
latter are more realistic since in the general theory it is unlikely that one will be able to Abelianize the constraints. 

We start from the classical expression (|49l) for the Hamiltonian constraint, 



H--^^ ^^^] -^+'^H^ + ar(^] =0. (89) 



(a; + 0)72 \ 8 I 2{x + a) 2Ev ' ' V^^*^ 
We then Abelianize the Hamiltonian (|49p . multiplying by ^^^"^"^ and grouping terms as 

H={^j^)'-l-^{i- + a)Aiy. (90) 

We wish to write the discretization in terms of classical quantities that are straightforward to represent in the 
quantum theory. Here one has to make choices, since there are infinitely many ways of discretizing a classical 
expression. In particular, we will notice that there exists, for this model, a way of discretizing the constraint in such 
a way that it remains first class (more precisely, Abelian) upon discretization. This is unusual, and we do not expect 
such a behavior in more general models. 

We now proceed to discretize this expression and to "holonomize" it, that is, to cast it in terms of quantities that 
are easily representable by holonomies. 



e 



{xm + a)'^e (xm-i + aye^ \ 1 



{{xrn + a) sin^ {pA^^„i) - {xm-i + a) siiY (pA^,,„„i)) 



{E^,r {El-,? J ^I'p' 

(91) 

expression that recovers (|90p in the limit e — > 0, p — > 0. From now on we will assume the spacing is uniform so we 
are dropping the m dependence of e. This expression is immediately Abelian since it can be written as the difference 
of two terms, one dependent on the variables at m and the other at m — 1. Therefore each term has automatically 
vanishing Poisson brackets with itself and with the other, 

H'n, = ^ {riXrr,,El,A^^„^) " ri^rn-l, E:^_„ A^,m-l)) , (92) 

with 



(l)''{Xrn,El,A^^rn) = ^^71^ \V ~ ^™ ~ a 9 9 (^"^ + o) siu^ ipA^^„i) . (93) 



At the horizon, the boundary condition is ^''(0,0, 0) = 0. The condition that the constraints vanish H^^ = is 
equivalent to (/)P{xm, E!^^, A^^m) = 0. 

When one has an Abelian set of constraints the uniform discretization quantization procedure is equivalent to the 
Dirac quantization procedure, as was shown in 8]. So for simplicity we will just follow the Dirac procedure. To 
implement the constraints as quantum operators as one does in the Dirac procedure, it is convenient to solve the 
constraints for the E^, ■ 



m 7 



= (94) 
1-^ + W^in^ (2p7if^,„^) 

and this relation can be immediately implemented as an operatorial relation and find the states that satisfy it. It 
should be noted that this relation can be implemented for other gauges as well in a straightforward manner. The 
states are given by, 

^[K^^^,T,a] =C(T,a)exp ±^3^ — J^J'^^'P.m] , (95) 

V *Planck ^ / 
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where C(t, a) is a function of the variables at the boundary r and a, which has to solve the constraint at the boundary 
(we will discuss this later). The functional / has the form, 



1 



-2F 



{xm + a)e 



4^2^2 I 



F sin(27piC^,„), 



472^2 1 




(96) 



with F{(f>, m) = Jq{1 — siv? t)^^/^dt the Jacobi Elliptic function of the first kind. Notice that the continuum limit 
of this expression for the state is immediate, i.e. the sum in m becomes an integral. 

We now need to impose the constraints on the boundary, in particular p'^ = —a/2 (in the limit N — > 00). Quantum 
mechanically p'^ = — J^ianck^/^''' ^^'^ therefore. 



C(t, a) — Co (a) exp 



2/'2 

^*-Planck 



(97) 



and Co (a) is an arbitrary function. This is analogous to the quantization that Kuchaf found where one had wave- 
functions that only depended on the mass. We have therefore completely solved the theory. 

Remarkably, the physical state we found is normalizable in a kinematic space of wavefunctions associated with one 
superselection sector, that is, the space of K'^s defined in the interval [0,7r/(7p)] (this defines an inner product in the 
bulk, the picture is easily completed in the boundary by considering functions Co (a) which are square integrable, as 
we will do later). In such space one can define operators associated to quantities that are not observables in the theory. 
This is therefore well suited to treat the problem of time in a relational way, as proposed by Page and Wootters [T^ . 
The main objection to the Page-Wootters construction was that one could not construct conditional probabilities 
based on states of the kinematical Hilbert space because the physical states were not contained in the kinematical 
Hilbert space and therefore the conditional probabilities were not well defined on physical states (see Kuchaf [2^ for 
a detailed discussion). Here, since the physical states exist in the kinematical space one does not face that problem. 
Although in the particular example we are considering there is not much point in defining a relational time, given its 
simple and static nature, it should be noted that the above result exhibits some level of robustness. For instance, it 
is true every time we are considering models where the variables are handled via the Bohr compactification. This is 
true, for instance, in cosmological models. It is also true if one couples the present model to a scalar field. It may also 
hold in more complicated models when one gauge fixes the diffeomorphisms. This is worthwhile further investigation. 

The above quantization has taken place in the "connection representation" . We would like to see that one obtains 
equivalent results in the loop representation. It is more convenient to write (j94p operatorially as Om"^ = where. 



472^2 



{x,n + a)e 



since this expression is straightforward to represent in the loop representation. 



1 - 



, + a 872/92 y (x^ + a)2e 



/;2 -72/4 

Pm ] -^Planck iTrf*,. \ 
^2^2 



u2 -y'^t 



Planck 



2/^^72^ 



,2,94 

Planck 



u2 -72* 



Planck 



2Mm7'^4lanck 



{Xm + a)2e2 {x,n + a)2e2 



*(^™ -4/9) 



1672/52 



(xr,i + a)2e2 (xra + a)^e^ 

= l'(Aim)- 



*(Ai,„ +4p) 



1672^2 



(99) 



The above equation is a recursion relation that implies that in the functions ^'(/im) the possible values of the /i's 
are /x = \x(r) — ±r + 4np with r € [0, 2p\ and n an integer. Solutions for different values of r are therefore not 
connected, and there is therefore a superselection. It is suggestive to compare this expression with the one obtained 
in loop quantum cosmology, where we see that for each point in the radial direction our wavefunctions have a similar 
recursion relation to those in loop quantum cosmology for the wavefunction of the universe. 

One can also see that the solution to this recursion relation can be obtained via the "loop transform" from the 
solution in the connection representation we found. Each solution in the connection representation corresponds to a 
given r-scctor of the superselection rule via. 



07 /-^/(P^) 

*r(/im) = — / dK^^{K^) exp (2/9 A" 
I" Jo 



(100) 
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The physical space of states is a Hilbert space with a natural inner product given by the norm in the variable 
a, and one must demand that the functions Co(a) be square integrable. The theory has two independent observables, 
the mass and its canonically conjugate variable. As mentioned above, the wavefunctions are functions of the mass. 
One can define observables with more geometric content, for instance the metric in a given gauge, as a function of 
the mass, using for example equations (|62II64[) . and similarly in Schwarzschild coordinates. The continuum limit is 
taken trivially, since the variables in a given coordinate system are uniquely defined in terms of the functions of the 
physical space and one approximate a continuum solution as accurately as one needs by reducing the stepsize. If one 
keeps the theory discrete, the relation is only approximate. It should be noted that the continuum limit is achieved 
in the limit e — s- and p — s- 0. If one adopts the point of view commonly used in loop quantum cosmology, that the 
quantum of distance should have a minimum value, then one would not expect to take the limits p and e going to 
zero, but to keep the parameters at a minimum value. In such a case one could expect to eliminate the singularity. 
This is plausible since then the triads would likely not go to zero. However, this would require a more careful analysis 
using coordinates that actually reach the singularity, which is not the case for the coordinates we have chosen so 
at the moment we cannot conclusively state what happens at the singularity. Also, it should be noted that in the 
most recent treatments of loop quantum cosmology the quantity that plays the same role as p is not taken to have 
a constant value, but to depend on the dynamical variables ^21j. If one followed a similar approach here, the results 
we derived would change. The results would remain valid in the asymptotic region far away from the horizon, since 
there p would tend to a constant. 

Another interesting aspect is that although we regularized the theory using a lattice, the level of ambiguity of the 
construction is limited, in part in order to have the constraints be Abelian. Although very limited, this example 
suggests that the use of a lattice regularization is not necessarily fraught with uncontrollable ambiguities. 

An aspect we have not emphasized is that in the calculations we have assumed that a > 0. Since a is the dynamical 
variable on which the wavefunctions end up depending, it is problematic to work with an operator with a continuum 
spectrum and positive eigenvalues only. A better way to handle this is to allow a to take all possible real values and 
then the gauge fixing condition we introduced should read E'^ ~ [x + |a|)^, and all subsequent equations involving 
X -\- a have to be modified accordingly (the conjugate momentum should be modified as well). In this treatment the 
mass is given by |a|/2 and is automatically positive. It is interesting to speculate what would happen if one wished to 
consider negative masses. What is clear at this stage is that the analysis would differ significantly with what we did 
in this paper, since in the Cauchy surface one would have to take into account the singularity and one cannot limit 
oneself to study the "exterior" as in the case of positive mass, where there is a well defined causal boundary at the 
horizon. The discussion of negative masses is better postponed till one can handle the interior problem and discuss 
the possibility of eliminating the singularity, since in the negative mass case the singularity has to be faced from the 
outset. 



IV. QUANTIZATION USING UNIFORM DISCRETIZATIONS 

Although we have succeeded in quantizing the model using the traditional Dirac quantization, we saw that in order 
to do this one had to use a property — the Abelianization of the Hamiltonian constraint — that is unlikely to hold 
in more general models. It is therefore worthwhile asking: what would have happened if we didn't Abelianize the 
constraint? In such case the traditional Dirac quantization would not have succeeded and we would have to resort 
to alternative proposals, like the uniform discretizations. Generically the discrete constraints that are obtained by 
discretizing a field theory with first class constraints are second class and become first class only in the continuum 
limit. We would like to develop a quantization strategy for such systems. An immediate answer to the problem is to 
deal with the second class constraints using the Dirac procedure. This is likely to be very onerous in cases of interest, 
and it does not take advantage of the fact that in the continuum limit the constraints become first class. 

Here is where our uniform discretization approach can help. In this approach the discrete theory has no constraints 
and nevertheless is capable of approximating the continuum theory in a controlled fashion. What we would advocate is 
to construct the discrete theory using the uniform discretization approach and then proceed to quantize the resulting 
discrete theory. In some cases one will be able to take the continuum limit in the quantum theory and this completes 
the quantization of the original continuum theory one started with satisfactorily. In some cases, as it is likely to be 
the case in the most interesting situations, it might occur that the continuum limit cannot be taken in the quantum 
theory. In such cases the approach we advocate is the following: what matters in a quantization procedure is to 
recover in the semiclassical limit the classical theory one started from plus corrections. We know the classical discrete 
theory we constructed approximates the continuum theory well. We therefore expect the quantum discrete theory to 
also approximate well the quantum continuum theory, even in cases where we cannot construct the latter exactly via 
this method, at least for certain states (it will obviously fail, for instance, for states that probe lengths smaller than 
the lattice spacing). We would therefore end with a quantum theory that approximates semiclassically the classical 
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theory we started with, plus corrections, i.e. the goal we were trying to achieve. 

Since these points of view imply a radical departure from the traditional Dirac method, it is worthwhile investigating 
how they perform in the face of concrete models where one can carry out the computations explicitly. In this section we 
would like to apply the uniform discretization approach to the treatment of spherically symmetric midi-superspaces. 
We will see that the approach is feasible but calculations can become quite involved if one does not choose to Abelianize 
the constraints. 

A. A brief summary of Uniform discretizations 

To recap briefly on previous discussions of uniform discretizations we would like to summarize its application to a 
field theory. We assume one is starting with a theory with variables q'^{x, t). One discretizes the underlying space-time 
manifold. The action can therefore be approximated by, 

S = ^L{n,n+l) (101) 

11 

L{n, n+l)^J2 ^^^iilkn], kkn+i]), (102) 

m 

where „ represents a discretization of q^{x,t) with m representing the array points in the spatial discretization 
and n the ones in the timelike direction. The notation [f?^ „] means the variable „ and some neighboring q'^'s 
determined by the scheme used to discretize spatial derivatives. For instance, if one uses a non-centered first order 
scheme one would have [9^ „] = (9m «' 9m+i «)■ ^ coordinate spatial separation between neighboring points 

in the lattice and A the time-like separation of neighboring points. Though we do not make it explicit, L depends 
on both e and A. We choose a first order approximation for the time derivatives since it simplifies constructing a 
canonical theory (for treatments with more than two time levels see (22|V One works out the canonical momenta for 
the discrete action, 

fc _ dL{n,n + l) 
Pm,n+i = -jTk • (103) 

One can also work out the momenta at n, which via the Lagrange equation, 

dL{n,n+l) dL{n—l,n 



yields. 



0, (104) 



pLn^-'-^t^. (105) 



These equations imply a canonical transformation between q,p at level n and level n + 1 via a type I generating 
function given by F{qn, qn+i) = —L{n,n + 1). It should be noted that we have not made a distinction between 
variables and Lagrange multipliers. Upon discretization, variables that in the continuum were Lagrange multipliers 
become evolution variables and they get determined by the evolution equations (this is a generic statement for 
Lagrange multipliers associated with diffeomorphism invariances, for many details concerning the Dirac treatment of 
discrete systems and particular examples see [23|). Equations that in the continuum were constraint equations now 
become evolution equations. The number of degrees of freedom is therefore larger than in the continuum theory. There 
will be different solutions in the discrete theory that correspond to different approximations or parameterizations of 
the same solution in the continuum theory. The determination of the Lagrange multipliers has proved problematic 
in previous examples we have studied. The resulting equations are usually polynomials of high order that can have 
complex solutions or develop discontinuities and branches. This may translate in that the discrete theory produces 
solutions that are not close to the constraint surface of the continuum theory throughout the whole evolution, and 
may depart significantly from it at certain points. This was usually considered a major obstacle to the use of these 
discrete theories. The uniform discretization approach bypasses this problem. 

The approach consists in replacing the evolution equations that we obtained for the dynamical variables of the 
problem by a set of evolution equations that in a sense we will make precise later on preserves the constraints of the 
continuum theory at a given level of approximation. The reason we can do this is that there is a large amount of 
freedom in how one discretizes a theory. In a sense we will be exploiting this freedom to our advantage. 
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To construct the uniform discretizations we note that since the evolution is given by a canonical transformation, 
generically one can find an infinitesimal generator for it. We will call this generator H. In terms of this generator, 
the discrete evolution of the dynamical variables can be written as, 

qLu+I = e^*'"^ {qln.n) ^ <C,n + {<lla.,nM + \{{ql,,,„MM + ' ' ' (106) 

We will assume a discretization has been chosen such that the form of H is given as H{q, p) = f{(j)i {q,p), . . . ,(t)p {q, p)) 
where (pi, . . . , 4>p are discretizations of the constraints of the continuum theory. To simplify the notation we have 
dropped the fact that the discretizations of the constraints will generically be functions of [g^j „], [pl^ „]. The function 

/ has the following properties: a) f{xi, . . . , xp) = ^=> Xa = 0, a = 1 . . . P and otherwise / > 0; b) ^'^^gj,'^''"'' — 0; 

c) det 7^ 0, Vx and d) f{(f>i{q,p), ■ ■ ■ , (t>p{q,p)) is defined for all {q^p) in the complete phase space. 

Discretizations of this sort can be constructed by suitable choices when one discretizes the action, as we proceeded. 
Further discussion on this is in reference (3]. However, we will not use in any way that the discretization was arrived 
to in this form so we refer the reader to the reference and just take the discretization as a given. 

A particularly simple example of such a function, that is useful in many systems is to just consider the sum of 
the squares of the discretized constraints. It is immediate that the evolution equations we wrote preserve exactly H. 
Therefore this means that if one chooses a small initial value for H one is guaranteeing that the sum square of the 
constraints is kept small and therefore one is not departing significantly from the constraint surface. This is also the 
case if one takes a more general / as well. 

Let us now show that this evolution does indeed have the correct continuum limit. For concreteness, and since 
this is general enough to accommodate the example we study in this paper, we consider a system that represents a 
discretization of a 1 + 1 dimensional field theory of the continuum with one first class constraint. We assume the spatial 
direction has a finite extension [a;o,a;i]. We have a lattice with N points such that \xi — xq] = Nt. The time-like 
direction starts at to and has spacing A. The field has R components q'^{x),i — 1 . . . R and canonical momenta p^{x). 
One has that. 





— q'^{xo + me,to + nA), 


(107) 


Pm,n 

e 


= p^xn + m€,to + nA), 


(108) 


X 


= lim xq + me. 


(109) 


t 


= lim Xq + toA, 


(110) 


q\x,t) 


= lim q^{xo + me,tQ + nA), 

e,A— >0 


(111) 


p\x,t) 


= lim p^{xq + me,tQ + nA). 


(112) 



We also assume a single (field-theoretic) constraint in the continuum that corresponds to N constraints 
4>a{[Pn,m], [<Zn,m],e) in the discrete theory. The generator is given by H = ^^Sa=i'^a/2- units where ft = c = 1 
the constraints have dimensions of length"^ and the constant k has dimensions of length so that the generator is 
dimensionless. We assume that the constraints have been appropriately rescaled such that. 



U2 

limH= — / dx<j)^[q{x),p{x)], (113) 

e^O 2 J 

where we emphasize that the constraint has a functional dependence on p, q. 

Since H is preserved upon evolution, we take H = (5^/2 with S a constant. We define Na = k<j>a{[q\, [p])/S which 
satisfies J2a=i — 1- We note that, 

^"'"^fc7^"'" = ^{9kn.II} + 0(A) =E{g™,„,0a([g],b])}A^a+O(A) (114) 

a— 1 

and identifying /S. — k5 and taking the limit e ^ (and therefore oo)^ (5-^0 one has that, 

q\x,t)^ [ dz{q\x,t),<j)[q{z,t),piz,t)]}N, (115) 
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where we have taken the hmit of vanishing time and spatial spacings and therefore the number of spatial points and 
of constraints to infinity, as one would expect in a field theory. In the discrete theory the canonical Poisson brackets 
are {qlnmPmn} = S^'''Sm,n with Kronecker deltas and in the continuum {q^ {x , t) , p' {z , t)} = 5^H{x — z). We have 
assumed that the model has one spatial dimension, as we consider in this paper, it is straightforward to generalize 
the construction to more spatial dimensions. Although we did not specify the field theory we were considering, it 
should be noted that in diffeomorphism invariant theories, to construct the generator one needs to integrate a density 
of weight one. Therefore one should choose appropriately densitized versions of the constraints. 

The above equations are valid if the constraints are first class in the continuum theory and the corresponding 
quantities in the discrete theory become first class constraints in the continuum limit. Gencrically, constraints that 
are first class in the continuum become second class quantities when discretized and will not necessarily become first 
class in the continuum limit of the discrete theory, which might also fail to exist altogether. To be precise, by "second 
class" in the discrete theory we mean that the vanishing of the quantities does not imply that their Poisson bracket 
vanishes. Notice that the equations specify the values of the Lagrange multipliers, they are given, up to a factor, by 
the values of the constraints. This therefore mandates that the constraints must be first class since only then one is 
free, in the continuum theory, to pick a particular set of Lagrange multipliers, which is what our method is doing. 
If the constraints are second class, the Lagrange multipliers are determined by the Dirac procedure and therefore 
the above procedure can be potentially inconsistent. There are two ways to proceed in this case. One of them is 
to use Dirac brackets in all the above expressions. Then the method still works as specified. Another possibility 
is to pretend the constraints are first class and then to study the continuum limit of the theory. If the continuum 
limit exists and in such limit the constraints become first class, then one may be able to use the discrete theory to 
construct the continuum limit and use it as a vehicle for quantization. Although it might appear that this latter point 
of view is less warranted than the first, in practice for cases of interest (as in general relativity in 3 + 1 dimensions) 
the use of the Dirac brackets will lead to unsolvablc equations and one may opt for the second avenue as a route for 
quantization. In fact, the latter point of view is likely to become the preferred one for cases of interest. The reason 
for this is that even in the case in which the continuum limit does not exist, one can use the constructed discrete 
theory as an approximation of the continuum one since the continuum constraints can be kept small. There might 
be other reasons why one does not want to consider the continuum limit. For instance, as occurs in loop quantum 
gravity, the kinematical space is not well suited for taking spatially continuum limits since one expects space to be 
quantized, as indicated by the quantization of the area, volume, etc. We will see how this operates in the example we 
discuss in this paper. 

B. Spherically symmetric vacuum gravity via uniform discretizations 

We would like to treat the model using uniform discretizations. Here one faces two possibilities. Since the discrete 
constraints are not first class, one can only ensure that the uniform discretizations will have the continuum limit if 
one treats them using Dirac brackets. However, one docs not expect that in more complicated models one will be able 
to take this path. We will therefore choose to treat the model using a technique more suited for more general models. 
The technique is to implement the uniform discretizations without using the Dirac brackets. Since one knows that in 
the continuum limit the constraints become first class, it is possible that in the continuum limit one will recover the 
desired theory. We will actually use some of the ambiguities in the discretization to make this outcome more likely. 

We will proceed to construct the quantity H which corresponds to the square of the Hamiltonian constraint. Here 
one faces a significant degree of ambiguity in choosing how to discretize the expression. To guide us in this we would 
like to require some conditions on the discretization that make it more likely that the continuum quantum limit will be 
achieved. The first condition is that when we take the continuum limit in the classical theory, the classical continuum 
constraint algebra should be reproduced. That is, the continuum limit of the Poisson bracket of the constraints of 
the discrete theory should reproduce the continuum classical constraint algebra. The second condition is that the 
spcctriim of the quantity H contain the zero eigenvalue in the combined classical and continuum limit. If these 
conditions are met, then one can be satisfied that one constructed a suitable quantum theory in the continuum limit. 

In achieving the above limit one expects two different types of problems. On the one hand, it could be that the 
constraints have been discretized in such a way that even classically they do not reproduce the constraint algebra in 
the continuum limit and this may imply that even classically H does not vanish in the continuum limit. In addition 
to this problem, one may have quantum anomalies. That is, even if one ensured that the algebra was reproduced 
classically in the continuum limit, upon quantization one may fail to reproduce the algebra in the continuum limit. 
The two types of problems can be characterized by the appearance of terms of order e" with e the spatial separation 
and n some power in the case of classical effects and the appearance of terms (^pianck/e)™ with m some power in the 
case of the quantum anomalies. For things to work out in the limit we should have that n > (m > due to how 
quantum anomalies arise). This suggests that the classical continuum limit should be taken ^pianck — > first and then 
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0. 

So what we are arguing is that due to the presence of quantum anomaUes one could face the situation where there is 
no well defined quantum continuum theory if m > 0, and nevertheless the discrete quantum theory is able to produce 
an acceptable continuum classical limit. This has important parallels with proposals by Klauder to extend the more 
traditional approaches to quantization [24] , and proposals on how to handle Thiemann's master constraint in the case 
in which it does not have zero eigenvalue [251] . 

We already have a candidate for the constraints that satisfies all the above conditions: the Abelianized version. 
What we would like to do, however, is to consider other versions of the constraint, where the above problems arise, 
and show that one can construct a suitable quantization. At this point one has to choose a discretization and a 
quantization of H that satisfies the two conditions we outlined above. There is no generic algorithmic procedure to 
construct such a discretization and quantization. In a generic situation one can imagine proposing parameterized 
discretizations and checking that the conditions are satisfied by choosing parameters. It is also possible to add terms 
that vanish in the continuum limit to the constraints to achieve the desired result. In the particular model we are 
studying we will construct a suitable discretization by starting from the Abelian version of the Hamiltonian constraint 
()90l) . and we multiply it times powers of the triad so one gets a non- Abelian constraint. In this case we will choose 
to multiply it times a fourth power of the triad since this simplifies quantizing it in the loop representation. The 
classical discrete Hamiltonian constraint we take is. 



HP = 



1 



3^ ,2 



{Xm + a) Sin^ {2p^K^^rn) ~ {Xm-l + a) S\V? {2pjK^^(^rn-l)) 



4-j,2p2 



The quantization we choose, in order for the operator to be Hermitian is, 
1 



HP = 



e'^{xm + a)^{xm-i + aY 



{El_,f {{xr,. + af + e') Y {Elf + af + e^) Y 

{xra + a) sin^ {2p-ik£^ - [xm-i + a) sin^ (2p-ikf^^_^ 



(116) 



872^2 



[xm + a) siv? (2pjk£j - (xm-i + a) sirY {2p^k'l_^^ 



8-f^p^ 



(117) 



From here one can construct the operator EI — (^HP^j , and we do not need to include powers of the determinant 

of the metric, since one has fixed spatial diffeomorphisms and therefore the spatial integral of the Hamiltonian as 
a scalar is well defined. The task now is to solve the eigenvalue problem and find the minimum eigenvalue of this 
positive definite quantum mechanical operator. This is a well defined problem in quantum mechanics, but it is far 
from trivial. Further studies on how to handle these types of problems are clearly needed. 

For the particular model at hand we can however determine an upper bound on the value of the minimum eigenvalue. 
To do this we use the states (|95p . which are normalized in the kinematical inner product 



— / dK^rn I da'i!*[K^^rn,T,a]^[K^rn,T,a]^l. 

^ ^ Jo Jo 



(118) 



We use this inner product to compute the expectation value of H for these states, and we get, 



(vI/]H]*) = 



(J ^Planck 



"^Planck 



8 

E 

n— 3 



n2n 

•^Planck 



(119) 



where the constants Ci are of order unity. Since the inner product has an integration over a the above expressions 
in the powers of a really correspond to expectation values, for instance the first term strictly speaking should be 
Cie^ {a~^) . The expectation value above satisfies the conditions we imposed for the continuum classical limit, i.e. it 
goes to zero as e ^ 0, ^pianck ~* 0- Moreover, in the continuum quantum limit (e ^ but ^pianck finite) it is divergent. 
Notice that this problem may arise even for the exact eigenvector corresponding to the minimum eigenvalue of H. 
That is, it is possible that one may not be able to use the discrete quantum theory to define a continuum quantum 
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theory as a limit. The minimum value of bound for the expectation value is achieved (assuming a large value of (a)) 
for e = \/^pianck(")- ^^^^ valuc (H) ^^pianck/i'*)^- ^^^^ for reference, for a Solar sized black hole, this amounts 
to 10-8°. 



V. CONCLUSIONS 



We have analyzed spherically symmetric quantum gravity in the connection and loop representations, outside the 
horizon. We were able to find a gauge fixing of the spatial diffeomorphisms that yields a theory with an Abelian Gauss 
law and a Hamiltonian constraint with a first class algebra with structure functions. A redefinition of the constraints 
yields an Abelian algebra of constraints. The latter can be quantized using the traditional Dirac procedure and yields 
a quantum theory in the connection and loop representation that coincides with Kuchaf 's treatment in terms of the 
traditional variables. We show one can discretize the theory to regularize the Hamiltonian during quantization in a 
reasonably unambiguous way and take the continuum limit, which is well defined. 

An interesting point is that in the gauge fixed model the physical quantum states of a given superselection sector 
are part of the kinematical Hilbert space. This allows, for instance, to implement in a well defined way a relational 
solution to the problem of time. This behavior is also expected in other models and needs to be further explored. 

Although we succeeded in implemented a traditional Dirac quantization for the model, this was done at a certain 
price. On one hand we gauge fixed the diffeomorphisms. On the other we Abelianized the Hamiltonian constraint 
through a rescaling. Gauge fixing the diffeomorphisms may be viewed as a legitimate approach to gain understanding 
into models, but ultimately is faced with the usual question of if the quantization obtained is equivalent to a quanti- 
zation performed in another gauge. The Abelianization of the Hamiltonian is a property that as far as we know, only 
holds in the gauge fixed model and appears unlikely to hold in other models generically (although an Abelianization 
without gauge fixing was achieved using the complex form of the Ashtekar variables in 

In view of the mentioned limitations of the traditional approach, we have also made a first exploration of the 
application of the uniform discretization approach to the system when we do not Abelianize the algebra of constraints. 
We have not done a complete analysis of such a model, since the calculations needed are significant. We note that 
it is plausible that one may produce a quantum theory that does not have a continuum limit, but has a classical 
continuum limit. This behavior may also be expected in more complicated models, where constraints generically 
cannot be Abelianized. Ambiguities in the discretization can be reduced by choosing things in such a way that the 
quantum theory approximates the continuum as much as possible. In this example, if one applied this criterium, one 
would recover the Abelian theory previously mentioned. 

Future developments will consist in studying the complete space-time using Kruskal-type coordinates and to see 
what effect one would have on the singularity due to the use of the loop representations. This should also be compared 
with mini-superspace treatments of the interior based on the isometry of the interior space-time with Kantowski-Sachs 
[26| . although the latter should be more limited since the isometry is limited to the interior region. We see in our 
treatment some initial hints of singularity removal since the discrete equations resemble those found in loop quantum 
cosmology, though further details need to be considered. 
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